SOLVING A GENERALIZED HERON PROBLEM 
BY MEANS OF CONVEX ANALYSIS 

BORIS S. MORDUKHOVICrfl NGUYEN MAU and JUAN SALINAS^ 

Abstract The classical Heron problem states: on a given straight line in the plane, find a point C 
such that the sum of the distances from C to the given points A and B is minimal. This problem can 
be solved using standard geometry or differential calculus. In the light of modern convex analysis, 
we are able to investigate more general versions of this problem. In this paper we propose and solve 
the following problem: on a given nonempty closed convex subset of IRf , find a point such that the 
sum of the distances from that point to n given nonempty closed convex subsets of 1R S is minimal. 

1 Problem Formulation 

Heron from Alexandria (10-75 AD) was "a Greek geometer and inventor whose writings 
preserved for posterity a knowledge of the mathematics and engineering of Babylonia, an- 
cient Egypt, and the Greco- Roman world" (from the Encyclopedia Britannica). One of the 
geometric problems he proposed in his Catroptica was as follows: find a point on a straight 
line in the plane such that the sum of the distances from it to two given points is minimal. 

Recall that a subset 0, of M s is called convex if Ax + (1 — X)y € f2 whenever x and y 
are in VL and < A < 1. Our idea now is to consider a much broader situation, where two 
given points in the classical Heron problem are replaced by finitely many closed and convex 
subsets Hi, i = 1, . . . , n, and the given line is replaced by a given closed and convex subset 
fl of M s . We are looking for a point on the set Q such that the sum of the distances from 
that point to fij, i = 1, . . . , n, is minimal. 

The distance from a point x to a nonempty set Q is understood in the conventional way 

d(x-,n) = M{\\x-y\\ {yen}, (1.1) 

where || • || is the Euclidean norm in ]R S . The new generalized Heron problem is formulated 
as follows: 

n 

minimize D(x) : = d(x; fl;) subject to x E Q, (1-2) 

i=l 

where all the sets f2 and f2j, i = 1, ...,n, are nonempty, closed, and convex; these are 
our standing assumptions in this paper. Thus (jl.2p is a constrained convex optimization 
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problem, and hence it is natural to use techniques of convex analysis and optimization to 
solve it. 



2 Elements of Convex Analysis 

In this section we review some basic concepts of convex analysis used in what follows. This 
material and much more can be found, e.g., in the books [21 [3j [5]. 

Let /: M s —> M := (— oo, oo] be an extended-real- valued function, which may be infinite 
at some points, and let 

dom f := {x e M s \ f(x) < oo} 
be its effective domain. The epigraph of / is a subset of M s x M defined by 

epi / := {(x,a) € M s+1 | x & dom / and a > f(x)}. 

The function / is closed if its epigraph is closed, and it is convex is its epigraph is a convex 
subset of M s+1 . It is easy to check that / is convex if and only if 

f(Xx + (1 - X)y) < A/0) + (1 - \)f{y) for all x,y € dom / and AG [0, 1]. 

Furthermore, a nonempty closed subset f2 of M s is convex if and only if the corresponding 
distance function f(x) = d(x;£l) is a convex function. Note that the distance function 
f(x) = d(x;£l) is Lipschitz continuous on M s with modulus one, i.e., 

I/O) -f(y)\ < \\x-y\\ for all x,yeM s . 

A typical example of an extended-real- valued function is the set indicator function 

f o if x g n, 

5(x;n):=\ ' (2.1) 

I oo otherwise. 

It follows immediately from the definitions that the set C M s is closed (resp. convex) if 
and only if the indicator function (|2.ip is closed (resp. convex). 

An element v € M s is called a subgradient of a convex function / : M s — > M at x G dom/ 
if it satisfies the inequality 

(v, x - x) < /O) - /O) for all x G M s , (2.2) 

where (•,•) stands for the usual scalar product in M s . The set of all the subgradients v 
in (I2.2p is called the subdifferential of / at x and is denoted by df(x). If / is convex and 
differentiable at x, then df(x) = {V/(x)}. 

A well-recognized technique in optimization is to reduce a constrained optimization 
problem to an unconstrained one using the indicator function of the constraint. Indeed, 
i G n is a minimizer of the constrained optimization problem: 

minimize f(x) subject to x € f2 (2-3) 
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if and only if it solves the unconstrained problem 

minimize f(x) + S(x; fi), x£M s . (2.4) 

By the definitions, for any convex function ip: ]R S — > M, 

x is a minimizer of ip if and only if £ d<p(x), (2-5) 

which is nonsmooth convex counterpart of the classical Fermat stationary rule. Applying 
(|2.5p to the constrained optimization problem (|2.3p via its unconstrained description (j2.4|) 
requires the usage of subdifferential calculus. The most fundamental calculus result of convex 
analysis is the following Moreau-Rockafellar theorem for representing the subdifferential of 
sums. 



Theorem 2.1 Let p^: ]R S — > M, i = 1, ...,m, be closed convex functions. Assume that 
there is a point x € n™ =1 domc^j at which all but (except possibly one) of the functions 
(pi, . . . ,<p m are continuous. Then we have the equality 

m m 
i=l i=l 

Given a convex set fi C M s and a point x6Q, the corresponding geometric counterpart 
of (|2.2p is the normal cone to at x defined by 



N{x;n) := {v € M s \ (v,x-x) < for all xeSl}. (2.6) 
It easily follows from the definitions that 

dS(x; Q) = N(x; Q) for every x € Q, (2.7) 



which allows us, in particular, to characterize minimizers of the constrained problem (J2J3J) 
in terms of the subdifferential (j2.2j) of / and the normal cone (|2.6p to 17 by applying 
Theorem 12.11 to the function <p(x) = f(x) + 5(x; £1) in (|2.5|) . 

Finally in this section, we present a useful formula for computing the subdifferential of 
the distance function (jl.ip via the unique Euclidean projection 

u(x-,n) -.= {x e n \ \\x-x\\ = d{x-,n)} (2.8) 

of x G M s on the closed and convex set C M s . 

Proposition 2.2 Let f2 7^ be a closed and convex of M s . Then 



' rx-n(x;fi) 



dd(x; Q) = < 



d(x; Q) 

N(x-,n)niB if xett 
where IB is the closed unit ball of M s . 



} if x <£ 



a 
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3 Optimal Solutions to the Generalized Heron Problem 

In this section we derive efficient characterizations of optimal solutions to the generalized 
Heron problem (|1.2p . which allow us to completely solve this problem in some important 
particular settings. 

First let us present general conditions that ensure the existence of optimal solutions to 
dT2]). 

Proposition 3.1 Assume that one of the sets and Q,i, i = 1, . . . ,n, is bounded. Then 
the generalized Heron problem f) 1 . 2 j) admits an optimal solution. 

Proof. Consider the optimal value 

7 := inf D(x) 

in (jl.2p and take a minimizing sequence {x^} C 0, with D(xk) — >■ 7 as k — > 00. If the 
constraint set O is bounded, then by the classical Bolzano- Weierstrass theorem the sequence 
{xk} contains a subsequence converging to some point x, which belongs to the set f2 due to 
it closedness. Since the function D{x) in (|1.2p is continuous, we have D(x) = 7, and thus 
X is an optimal solution to (jl . 2[) . 

It remains to consider the case when one of sets Jlj, say fii, is bounded. In this case we 
have for the above sequence {x^} when k is sufficiently large that 

d(x fc ;fii) < D{x k ) < 7 + 1, 

and thus there exists wj, G with — Wk\\ < 7 + 1 for such indexes k. Then 

\\xk\\ < \\wk\\ + 7 + ! 5 

which shows that the sequence {x k } is bounded. The existence of optimal solutions follows 
in this case from the arguments above. A 

To characterize in what follows optimal solutions to the generalized Heron problem (|l,2p , 
for any nonzero vectors u, v £ M s define the quantity 

cos(v,u) := ■ (3.1) 

\\v\\ ■ \\u\\ 

We say that $7 has a tangent space at x if there exists a subspace L = L(x) ^ {0} such that 

N(x;Q) = L L := {v G M s \ (v,u) = whenever u G L). (3.2) 

The following theorem gives necessary and sufficient conditions for optimal solutions to 
(jl.2p in terms of projections (I2.8P on fij incorporated into quantities (13. ip . This theorem 
and its consequences are also important in verifying the validity of numerical results in the 
Section 3. 
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Theorem 3.2 Consider problem (jl.2p in which 

fii n = /or aZZ i = 1, . . . , n. (3.3) 

Given x € define the vectors 

x — n(x; Oi) . . 

' V O, i = l,...,n, 3.4 
ct(x; SZjj 

XTien x £ Q is an optimal solution to the generalized Heron problem (jl.2p i/ and on/y i/ 

- J^ai(x) € JV(x;fi). (3.5) 

i=l 

Suppose in addition that the constraint set £2 has a tangent space L at x. Then 113. 5\) is 
equivalent to 

n 

cos (oi(rr), n) = whenever u € L \ {0}. (3-6) 

i=l 

Proof. Fix an optimal solution x to problem (|1 .2[) and equivalently describe it as an optimal 
solution to the following unconstrained optimization problem: 

minimize D{x) + 6(x; D,), x € M s . (3-7) 

Applying the generalized Fermat rule (|2.5p to (|3.7p . we characterize x by 

n 

O€0(5^d(-;fi<) + 5(-;fi))(z)- (3.8) 
i=i 

Since all of the functions d(-; i = 1, . . . ,ra, are convex and continuous, we employ the 
subdifferential sum rule of Theorem 12.11 to (|3.8p and arrive at 

n 

o g <9(l> + 5(; n)) (x) = dd (x; n<) + iv(x; o) 

(3.9) 

= J^ai(x) + N(x;Q), 
i=l 

where the second representation in (|3.9p is due to (|2.7p and the subdifferential description of 
Proposition 12.21 with ai(s) defined in (|3.4p . It is obvious that (|3.9p and (|3.5p are equivalent. 

Suppose in addition that the constraint set O has a tangent space L at x. Then the 
inclusion (|3.5p is equivalent to 

<G J^fli(x) 
i=l 

which in turn can be written in the form 



a»(x), u ) = for all u € L. 



i=l 
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Taking into account that ||aj(x)|| = 1 for alH = 1, . . . ,n by (|3.4p and assumption (|3.3p . the 
latter equality is equivalent to 

y {a ;W' V ] „ = for all u £ L\ {0}, 
jri \\<h(?)\\ ■ \\u\\ 

which gives (13, 6p due to the notation (|3.1|) and thus completes the proof of the theorem. A 

To further specify the characterization of Theorem 13.21 recall that a set A of JR S is an 
affine subspace if there is a vector a G A and a (linear) subspace L such that A = a + L. In 
this case we say that A is parallel to L. Note that the subspace L parallel to A is uniquely 
defined by L = A — A = {x — y \ x £ A, y £ A} and that A = b + L for any vector 6 G A 

Corollary 3.3 Let f2 6e an affine subspace parallel to a subspace L, and let assumption 
(|3.3p of Theorem \3.2\ be satisfied. Then x G is a solution to the generalized Heron problem 
(jl.2p i/ and onZy «/ condition (|3.6p holds. 

Proof. To apply Theorem 13.21 it remains to check that L is a tangent space of f2 at x 
in the setting of this corollary. Indeed, we have = x + L, since is an affine subspace 
parallel to L. Fix any v G iV(x; J7) and get by (|2.6p that (v, x — x) < whenever i£!l and 
hence {v, u) < for all it G L. Since L is a subspace, the latter implies that (t>, u) = for 
all u £ L, and thus -/V(x; f2) C L . The opposite inclusion is trivial, which gives (13.2p and 
completes the proof of the corollary. A 

The underlying characterization (j3.6|) can be easily checked when the subspace L in 
Theorem 13.21 is given as a span of fixed generating vectors. 

Corollary 3.4 Let L = span{ui, . . . , u m } with Uj ^ 0, % = 1, ...,m, in the setting of 
Theorem 13. 21 Then x £ Vt is an optimal solution to the generalized Heron problem (|1.2p if 
and only if 

n 

cos (ai(x), uj) = for all j = 1, . . . , m. (3.10) 

i=l 

Proof. We show that (j3.6|) is equivalent to (|3.10p in the setting under consideration. Since 
(|3.6p obviously implies (|3.10p . it remains to justify the opposite implication. Denote 

n 

a := yaj(x) 
i=i 

and observe that (I3.10P yields the condition 

(a, Uj) = for all j = 1, . . . m, (3-H) 
= 1 for all i = 1, . . . , n. Taking now any vector 

with some Xj G M n 



since Uj ^ for all j = 1, . . . , m and ||aj 
u G L \ {0}, we represent it in the form 

m 

3=1 
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and get from (|3.1ip the equalities 



n 

(a,u) = ^2\j{a,Uj) = 0. 

This justifies (|3.6p and completes the proof of the corollary. A 

Let us further examine in more detail the case of two sets f^i and 0,2 in f)1.2|) with 
the normal cone to the constraint set VL being a straight line generated by a given vector. 
This is a direct extension of the classical Heron problem to the setting when two points are 
replaced by closed and convex sets, and the constraint line is replaced by a closed convex 
set Q with the property above. The next theorem gives a complete and verifiable solution 
to the new problem. 

Theorem 3.5 Let Q\ and ^2 be subsets of M s as s > 1 with Q n fij = for i = 1,2 in 
(jl.2p . Suppose also that there is a vector a^O such that N(x; fi) = spanja}. The following 
assertions hold, where a, := ai{x) are defined in (|3,4p : 

(i) IfxEQ is an optimal solution to (jl.2p . then 

either a\ + a<i = or cos (ai, a) = cos (02, a). (3-12) 

(ii) Conversely, if s = 2 and 

either a\ + 02 = or [017^02 and cos(ai, a) = cos(a2, a)] , (3.13) 
then x £ £1 is an optimal solution to the generalized Heron problem (jl.2p . 

Proof. It follows from the above (see the proof of Theorem I3.2[) that x € is an optimal 
solution to (jl.2p if and only if —a\ — G iV(x; fi). By the assumed structure of the normal 
cone to O the latter is equivalent to the alternative: 

either a\ + ai = or a\ + a<i = Xa with some A 7= 0. (3-14) 

To justify (i), let us show that the second equality in (|3.14p implies the corresponding 
one in (|3.12|) . Indeed, we have ||ai|| = ||oi|| = 1, and thus (|3.14p implies that 

A 2 ||a|| 2 = ||oi + a 2 || 2 = ||ai|| 2 + ||a 2 || 2 + 2(01,02) = 2 + 2(ai,a 2 ). 

The latter yields in turn that 

(ai, Xa) = (Xa — 02, Aa) 

= A 2 | |a| | 2 — A(a2, a) 

= 2 + 2(01,02) - A(a 2 , a) 

= 2(a 2 ,a 2 ) + 2(ai,a 2 ) - A(a 2 ,a) 

= 2(a 2 + 01,02) - A(a 2 ,a) 

= 2(Aa, 02) — A(o2, a) = (02, Aa), 
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which ensures that (a±,a) = (a 2 ,a) as A / 0. This gives us the equality cos(ai,a) = 
cos(a2,a) due to ||ai|| = ||a2|| = 1 and a ^ 0. Hence we arrive at f|3. 12|) . 

To justify (ii), we need to prove that the relationships in f|3. 13j) imply the fulfillment of 

— a\ — a 2 € N(x; Q) = span{a}. (3.15) 

If — ai — a2 = 0, then f)3. 15|) is obviously satisfied. Consider the alternative in (|3. 13|) when 
o-i 7^ 0-2 and cos(ai,a) = cos(a2,a). Since we are in 1R 2 , represent a\ = (x±,yi), a 2 = 
(x 2 ,y 2 ), and a = (x,y) with two real coordinates. Then by (|3.ip the equality cos(ai,a) = 
cos(a2,a) can be written as 

xix + yiy = x 2 x + y 2 y, i.e., (x± - x 2 )x = (y 2 - yi)y. (3.16) 

Since a 7^ 0, assume without loss of generality that y 7^ 0. By 

11 i|2 11 ||2 , . 2 1 2 2,2 

||oi|| =\\a2\\ x 1 +y 1 = x 2 + y 2 

we have the equality [x\ — x 2 ){x\ + x 2 ) = (y 2 — yi)(y 2 + yi), which implies by (|3.16p that 

y(xi - x 2 )(xi + x 2 ) = x{xi - x 2 )(y 2 + yi). (3.17) 

Note that x\ 7^ x 2 , since otherwise we have from (|3.16p that y\ = y 2 , which contradicts the 
condition a\ 7^ a 2 in (|3.13p . Dividing both sides of (|3.17j) by x\ — x 2 , we get 

y(xi +x 2 ) = x(y 2 +1/1), 

which implies in turn that 

y(ai + a 2 ) = y(xi + x 2 ,yi + y 2 ) = (x(yi + y 2 ), y(yi + y 2 )) = (yi + y 2 )a. 
In this way we arrive at the representation 

yi + V2 
a\ + a 2 = a 

y 

showing that inclusion (|3.15p is satisfied. This ensures the optimality of x in (|1.2p and thus 
completes the proof of the theorem. A 

Finally in this section, we present two examples illustrating the application of Theo- 
rem 13.21 and Corollary 13.41 respectively, to solving the corresponding the generalized and 
classical Heron problems. 

Example 3.6 Consider problem (jl.2p where n = 2, the sets fii and £l 2 are two point A and 
B in the plane, and the constraint O is a disk that does not contain A and B. Condition 
(|3.5|) from Theorem 13.21 characterizes a solution M £ SI to this generalized Heron problem 
as follows. In the first case the line segment AB intersects the disk; then the intersection 
is a optimal solution. In this case the problem may actually have infinitely many solutions. 
Otherwise, there is a unique point M on the circle such that a normal vector it to Q at M 
is the angle bisector of angle AMB, and that is the only optimal solution to the generalized 
Heron problem under consideration; see Figure 1. 



S 



Figure 1: Generalized Heron Problem for Two Points with Disk Constraint. 



Example 3.7 Consider problem (jl,2p . where f2j = {^4i}, i = 1, . . . ,n, are n points in the 
plane, and where O = £ C M 2 is a straight line that does not contain these points. Then, 
by Corollary 13.41 of Theorem 13. 2\ a point M 6 C is a solution to this generalized Heron 
problem if and only if 

cos(MA[, ~ct) H h cos(MA n , "S") = 0, 

where a is a direction vector of C. Note that the latter equation completely characterizes 
the solution of the classical Heron problem in the plane in both cases when A\ and Ai are 
on the same side and different sides of C; see Figure 2. 




4 Numerical Algorithm and Its Implementation 

In this section we present and justify an iterative algorithm to solve the generalized Heron 
problem (]1.2p numerically and illustrate its implementations by using MATLAB in two 
important settings with disk and cube constraints. Here is the main algorithm. 
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Theorem 4.1 Let Q and Qi, i = 1,... ,n, be nonempty closed convex subsets of M s such 
that at least one of them is bounded. Picking a sequence {a^} of positive numbers and a 
starting point x\ E £1, consider the iterative algorithm: 

n 

Xk+i = n(x fc - a k ^2v ik ;n\ A; = 1,2,..., (4.1) 

i=l 

where the vectors vu~ in (|4.ip are constructed by 

Vik '■= 37 tH with ^ik '■= n(x fc ;Oj) if x k £ fi* (4.2) 

d{Xk;Ui) 

and Vik := otherwise. Assume that the given sequence {ct k } in (|4.ip satisfies the conditions 

oo oo 

^^afe = oo and ^^af. <oo. (4.3) 
fc=i fc=i 

Then the iterative sequence {xk} in (|4.2p converges to an optimal solution of the generalized 
Heron problem (jl.2p and the value sequence 

V k :=min {D( Xj )\ j = 1, . . . , k) (4.4) 

converges to the optimal value V in this problem. 

Proof. Observe first of all that algorithm (|4.1|) is well posed, since the projection to a convex 
set used in (I4.2p is uniquely defined. Furthermore, all the iterates {xk} in (14. ip are feasible; 
see the proof of Proposition 13. 11 This algorithm and its convergence under conditions (j4.3|) 
are based on the subgradient method for convex functions in the so-called "square summable 
but not summable case" (see, e.g., [1]), the subdifferential sum rule of Theorem 12. 1\ and 
the subdifferential formula for the distance function given in Proposition 12.21 The reader 
can compare this algorithm and its justifications with the related developments in [3] for 
the numerical solution of the (unconstrained) generalized Fermat-Torricelli problem. A 

Let us illustrate the implementation of the above algorithm and the corresponding cal- 
culations to compute numerically optimal solutions in the following two characteristic ex- 
amples. 

Example 4.2 Consider the generalized Heron problem (11.2p for pairwise disjoint squares 
of right position in M 2 (i.e., such that the sides of each square are parallel to the x-axis or 
the y-axis) subject to a given disk constraint. Let q = (aj,6j) and r^, % = 1, . . . , n, be the 
centers and the short radii of the squares under consideration. The vertices of the ith square 
are denoted by qu = (a, + n, bi + n), q 2 i = (a, - n, h + r^), q 3i = (a; - r i; 6j - r^), q^ = 
(aj + rj, bi — rj). Let r and p = (^,77), be the radius and the center of the constraint. Then 
the subgradient algorithm (|4.ip is written in this case as 

n 

i=i 
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where the projection P(x,y) := H((x,y);£l) is calculated by 

r(x — v) 



P(x, y) = (w x +u, Wy+rj) with w x 



^{x - v) 2 + (y -rf) 
The quantities in the above algorithm are computed by 

if \x lk — Oi| < Ti and \x 2k - h\ < r, 



r(y - rf) 

_ and w y = — ^ 

y/(x - u) 2 + (y — rf) 2 



Xk - qu 



Vik = < 



\\Xk 


- fell 


Xk 


- fe 


\\x k 


- fell 


Xk 


- fe 


II** 


- fell 




- fe 


\\xk 


- fell 



(0,1) 



if x xk —ai>n and x 2 k - h > n, 

if xik —ai< -Ti and x 2 k - h > n, 

if x^ - a.i < ~n and x 2 k - h < -r, 

if x^ —ai>n and x 2 fc - h < -n, 
if |xi fc - a,-! < r i and x 2k - 6 l > r i} 



(0, -1) if \xik - cii\ < n and x 2k — < — r,-, 



(1,0) 



if [rcift - Oi| > n and \x 2k ~h\< n 



(—1,0) if \xi k - Oi| < -r» and |x 2 fc - b* | < n 
for all i = 1, . . . , n and k = 1,2,... with the corresponding quantities defined by (14. 4j) . 




-4 -2 2 4 6 



MATLAB RESULT 



k 


Xk 


v k 


1 


(-3,5.5) 


30.99674 


10 


(-1.95277,2.92608) 


26.14035 


100 


(-2.02866,2.85698) 


26.13429 


1000 


(-2.03861,2.84860) 


26.13419 


10,000 


(-2.03992,2.84750) 


26.13419 


100,000 


(-2.04010,2.84736) 


26.13419 


200,000 


(-2.04011,2.84735) 


26.13419 


400,000 


(-2.04012,2.84734) 


26.13419 


600,000 


(-2.04012,2.84734) 


26.13419 



Figure 3: Generalized Heron Problem for Squares with Disk Constraint. 

For the implementation of this algorithm we develop a MATLAB program. The following 
calculations are done and presented below (see Figure 3 and the corresponding table) for 
the disk constraint Q with center (—3,4) and radius 1.5, for the squares fit with the same 
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short radius r = 1 and centers (—7,1), (—5,-8), (4,7), and (5,1), for the starting point 
x\ = (—3,5.5) € il, and for the sequence of = 1/k in (|4.ip satisfying conditions (|4.3p . 
The optimal solution and optimal value computed up to five significant digits are x = 
(-2.04012,2.84734) and V = 26.13419. 

The next example concerns the generalized Heron problem for cubes with ball constraints 
in M 3 . 

Example 4.3 Consider the generalized Heron problem (jl.2p for pairwise disjoint cubes of 
right position in if? 3 subject to a ball constraint. In this case the subgradient algorithm 
jUD is 

n 
i=l 

where the projection H((x,y, z);£l) and quantities Vik are computed similarly to Exam- 
ple ED 




MATLAB RESULT 


k 


Xk 


v k 


1 


(2,2,0) 


27.35281 


1,000 


(-0.68209,0.25502,0.69986) 


24.74138 


1,000,000 


(-0.77641,0.31416,0.74508) 


24.73757 


2,000,000 


(-0.77729,0.31480,0.74561) 


24.73757 


3,000,000 


(-0.77769,0.31509,0.74584) 


24.73757 


3,500,000 


(-0.77782,0.31518,0.74592) 


24.73757 


4,000,000 


(-0.77792,0.31526,0.74598) 


24.73757 


4,500,000 


(-0.77801,0.31532,0.74604) 


24.73757 


5,000,000 


(-0.77808,0.31538,0.74608) 


24.73757 



Figure 4: Generalized Heron Problem for Cubes with Ball Constraint. 

For the implementation of this algorithm we develop a MATLAB program. The Figure 4 
and the corresponding figure present the calculation results for the ball constraint f2 with 
center (0,2,0) and radius 2, the cubes J7j with centers (0,-4,0), (6,2,-3), (—3,-4,2), 
(—5, 4, 4), and (—1, 8, 1) with the same short radius r = 1, the starting point x\ = (2, 2, 0), 
and the sequence of = 1/k in (|4.ip satisfying (|4.3p . The optimal solution and optimal 
value computed up to five significant digits are x = (—0.77808,0.31538,0.74608) and V = 
24.73756. 
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